home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libconv / TRY / simple_diir1d.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  332 b   |  21 lines
  1.     subroutine simple_diir1d( nf, ng, f, g, h)
  2.     double precision f(0:(nf-1)), g(0:(ng-1)), h(0:(nf+ng-2))
  3.     integer nf, ng
  4.  
  5.     integer i, j, k
  6.  
  7.     do k = 0, nf-1
  8.         h(k) = f(k)
  9.     end do
  10.  
  11.     do j = 0, nf-1
  12.         k = max( j - (ng-1), 0)
  13.         do i = k, j-1
  14.         h(j) = h(j) - h(i) * g(j-i)
  15.         end do
  16.         h(j) = h(j) / g(0)
  17.     end do
  18.  
  19.     return
  20.     end
  21.